Phase diagram for rotating compact stars with two high density phases 
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For the classification of rotating compact stars with two high density phases a phase diagram in the 
angular velocity (Q) - baryon number (N) plane is investigated. The dividing line N CI it(Q.) between 
configurations with one and two phases is correlated to a local maximum of the moment of inertia 
and can thus be subject to experimental verification by observation of the rotational behavior of 
accreting compact stars. Another characteristic line, which also can be measured is the transition 
line to black holes that of the maximum mass configurations. The positions and the shape of these 
lines are sensitive to changes in the equation of state (EoS) of stellar matter. A comparison of the 
regional structure of phase diagrams is performed for polytropic and relativistic mean field type 
EoS and correlations between the topology of the transition lines and the properties of two-phase 
EoS are obtained. Different scenarios of compact star evolution are discussed as trajectories in the 
phase diagram. It is shown that a population gap in the Q — N plane for accreting compact stars 
would signal a high density phase transition and could be observed in the distribution of so called 
Z sources of quasi periodic oscillations in low-mass X-ray binaries. 
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At present, the existence of exotic phases of matter at 
high densities is under experimental investigation in ul- 
trarelativistic heavy-ion collisions jjj the most prominent 
being the deconfined phase of QCD |Q. While the diag- 
nostics of a phase transition in experiments with heavy- 
ion beams faces the problems of strong noncquilibrium 
and finite size, the dense matter in a compact star forms 
a macroscopic system in thermal and chemical equilib- 
rium for which effects signalling a phase transition shall 
be most pronounced. 

Signals of high density phase transitions have been sug- 
gested in the form of characteristic changes of observables 
such as the surface temperature J3|,|| , brightness Q , pulse 
timing Q and rotational mode instabilities during the 
evolution of the compact object. In particular the pulse 
timing signal has attracted much interest since it is due 
to changes in the kinematics of rotation. Thus it could 
be used not only to detect the occurrence, but also to 
determine the size of the high density matter core from 
the magnitude of the braking index deviation from the 
magnetic dipole value Besides of the isolated pul- 
sars, one can consider also the accreting compact stars in 
low-mass X-ray binaries (LMXBs) as objects from which 
we can expect signals of a high density phase transition 
in their interior . The observation of quasiperiodic 
brightness oscillations (QPOs) for some LMXBs has 
lead to very stringent constraints for masses and radii 
jl2| which according to jl3) could even favour strange 
quark matter interiors over hadronic ones for these ob- 



jects. Due to the mass accretion flow these systems are 
candidates for the formation of the most massive com- 
pact stars from which we expect to observe signals of the 
transition to either quark core stars, to a third family of 
stars jF|| or to black holes. 

In this work we introduce a classification of rotating 
compact stars in the plane of their angular frequency O 
and mass (baryon number N) which we will call phase 
diagram. In this diagram, configurations with high den- 
sity matter cores are separated from conventional ones 
by a critical phase transition line. 

Since the phase diagram of rotating compact objects 
seems to be a more general approach for investigations 
of phase transition effects in the interior of the star we 
assume that the deconfinement transition could be a par- 
ticular case besides of other possibilities for phase transi- 
tions like pion or kaon condensation as discussed, e.g. in 
p5| , ^6| . Therefore, our aim is to suggest it as a heuristic 
tool for obtaining constraints on the EoS at high densities 
from the rotational behaviour of compact stars. We will 
provide criteria under which a particular astrophysical 
scenario with spin evolution could be qualified to signal 
of high density phase transition. 

The true EoS that describes the interior of compact 
stars is largely unknown. This results from the inability 
to verify experimentally the different theories that de- 
scribe the strong interactions between baryons and the 
many-body theories of dense matter at densities larger 
than about twice the nuclear density |l7|]. 
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Since our focus is on the elucidation of qualitative fea- 
tures of signals from the high density phase transition in 
the pulsar timing we will use a generic form of an equa- 
tion of state (EoS) with such a transition. We use the 
polytropic type equation of state Jig] in the form 
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where n is the baryon number density, Pj - the pres- 
sure, and 6i- energy density for both the low (i = L) 
and the high (i = H) density phases, respectively. = 
Ki(no) is the value of the incompressibility |25| Ki(n) = 
9 dP/dn at the saturation density no = 0.17//m 3 , I\ = 
d\n{Pi)/dhxei the adiabatic index and m is the nucleon 
mass. The phase transition between the lower and higher 
density phases is made by the Maxwell construction ]l9| ] 
and compared to a relativistic mean field model consist- 
ing of a linear Walecka plus dynamical quark model EoS 
H|,||| with a Gibbs construction |^J|, see Fig{| 
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FIG. 1. Upper panel: Two polytropic type EoS models: 
SH-S (dashed line) and H-SH (dotted line) are compared to a 
relativistic mean field (RMF) EoS, which has a hard-soft-hard 
(H-S-H) phase transition between a linear Walecka model 
for the low density hadronic phase and a dynamical quark 
model for the high density quark matter phase with a soft 
mixed phase j^]. SH-S (Tl = 2.5, Th = 1.9) represents a 
semihard-soft phase transition; H-SH (r_L = 2.6, Th = 2.2) 
hard-semihard one. Lower panel: Incompressibilities for the 
EoS model of the upper panel (same line styles) in units of 
the value RMF in (3- equilibrium, K° RMF = 1125 MeV. The 
shaded region separates the hard (H) and soft (S) parts of 
EoS and defines the semihardness (SH). 

The quark matter part of this EoS is obtained from a 
dynamical confining approach pl[ | in the generalization 
to three flavors p2], where the strange flavor remains 



confined at the deconfinement transition for the light and 
appears only at densities for which stars are close to the 
gravitational instability. The difference to most of the 
models for quark deconfinement in neutron star matter 
is that the ambiguity in the choice of the bag constant for 
the quark matter phase can be removed by a derivation 
of this quantity [ p3|j24| within the dynamical confining 
approach 
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With the EoS models discussed above we have per- 
formed calculations of rotating compact star configura- 
tions assuming stationary, rigid rotation. For our treat- 
ment of rotation within general relativity we employ a 
perturbation expansion following Refs. |26). For a de- 
tailed discussion of the method and its application we 
refer to ||. The results of our calculations of rotating 
compact star configurations can be classified in the plane 
of angular velocity Q and baryon number N which we call 
phase diagram. 
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FIG. 2. Phase diagrams for the RMF EoS (upper panel) 
and two polytropic model EoS (lower panels). 

In Fig. H we display the phase diagrams for the ro- 
tating star configurations, which correspond to the three 
model EoS of Fig. 0. These phase diagrams have four 
regions: (i) the region above the maximum frequency 
O > Q.k(N) where no stationary rotating configurations 
are found, (ii) the region of black holes N > A^ max (0), 
and the region of stable compact stars which is subdi- 
vided by the critical line iV cr it(fi) into (iii) the region of 
hybrid stars for N > iV cr i t (fi) where configurations con- 
tain a core with a second, high density phase and (iv) 
the region of mono-phase stars without such a core. 

From the comparison of the regional structure of these 
three different phase diagrams in Fig. |^ with the corre- 
sponding EoS in Fig. [I] we arrive at the main result of 
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this paper that there are the following correlations be- 
tween the topology of the lines N max (Cl) and iV cr it(fi) 
and the properties of two-phase EoS: 

- The hardness of the high density EoS determines 
the maximum mass of the star, which is given by 
the line N max (Cl). Therefore A r max (0) is propor- 
tional to the parameter Kh(iih), where nn is the 
density of the transition to high density phase. 

- The onset of the phase transition line iVcrit(O) de- 
pends on the density n# and Kl(jil) where rih 
is the density of the transition to the low density 
phase. 

- The curvature of the lines N max (Cl) and N cl -i t (ft) is 
proportional to the compressibility of the high and 
low density phases, respectively. 

Therefore, a verification of the existence of the critical 
lines iV cr it(r2) and -/V max (f2) by observation of the rota- 
tional behavior of compact objects would constrain the 
parameters of the EoS for neutron star matter. 

The determination of N max (Cl), the border between 
compact stars and black holes, is a traditional issue which 
has recently gained new impetus due to the interpreta- 
tion of recent LMXB data. A new aspect characterizing 
the configurations with a phase transition is the criti- 
cal line N CI it(Q) which can be measurable by changes in 
the rotational dynamics since it is correlated with a lo- 
cal maximum of the moment of inertia I(N, CI), the key 
quantity governing the rotational evolution via 



CI = 



K(N, CI) 



(1) 



i(N,n) + n(di(N,n)/dci) N ' 

In Fig.|| we show for the example of the RMF EoS with 
a deconfinement phase transition from hadronic to quark 
matter how the structure of the lines of constant moment 
of inertia correlated with the critical lines in the phase 
diagram. 

In Eq. ([!]) K — K int + K cxt is the net torque acting on 
the star due to internal and external forces. The internal 
torque is given by K int (N, CI) = -CLN (dI(N, Cl)/dN) n , 
the external one can be subdivided into an accretion and 
a radiation term K cxt = K acc + -fsT ra d- The first one is 
due to all processes which change the baryon number, 
K acc — N dJ/dN and the second one contains all pro- 
cesses which do not. For the example of magnetic dipole 
and/or gravitational wave radiation it can be described 
by a power law if ra d = /3Cl n , see [E7 28 1. 

To prove that the appearance or disappearance of a 
high density phase during the rotational evolution of the 
star could entail observational consequences for the an- 
gular velocity we consider three main representatives dif- 
ferent classes of trajectories which could cross the critical 
line on phase diagram. These classes of tracks are: (a) 
spindown of isolated (non-accreting, N = 0) pulsars due 



to magnetic dipole radiation [|6|j8j , (b) spin up in accreting 
systems with weak magnetic field |9|,|l0| (N ~ const, ver- 
tical tracks) and (c) accretion either with strong magnetic 
field p0[ o r for accreting binaries emitting gravitational 
waves p9| , for which CI ~ const (horizontal tracks). 
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FIG. 3. Phase diagram for configurations of rotating com- 
pact objects in the plane of angular velocity CI and mass 
(baryon number TV). Contour lines show the values of the 
moment of inertia in MqW. The line iVcrit(fi) which sep- 
arates hadronic from quark core stars corresponds the set of 
configurations with a central density equal to the critical den- 
sity for the occurrence of a pure quark matter phase. 

In the case of (a) the spindown K — K ra d or (b) spinup 
K = K acc evolutions (in both cases K- lnt « K ext ) the 
objects can undergo a phase transition if the baryon num- 
ber lies within the interval N CI i t (Cl = 0) < N < N c , 
where N c is the end piont of the critical line N cr i t (Cl). 

If the core of compact star is soft enough (as in case 
(SH-S)) the critical line N CI i t (Ct) crosses N max (Cl) at some 
Cl c (N c ) < CIk(N c ). This means, that massive mono- 
phase configurations with total baryon number N > 7V C 
rotating with angular velocities in the interval Cl c < CI < 
CIk should encounter a transition to a black hole during 
the spin down evolution. 

In cases, when the core EoS is harder (H-S-H and SH- 
H), the region of hybrid stars is a band which separates 
mono-phase configurations from black holes, see upper 
two panels of Fig. ^. 

As it has been shown in || for the vertical tracks (a) 
and (b) in the phase diagram, the braking index n(Cl) 
changes its value from n(Cl) > 3 in the region (iii) to 
n(Cl) < 3 in (iv). This is the braking index signal for a 
deconfinement transition introduced in Ref. ||. 

The third evolutionary track is accretion with strong 
magnetic fields [ [l0[ and/or gravitational wave emission 
(horizontal tracks) p9|.For this case the CI first decreases 
as long as the moment of inertia monotonously increases 
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with TV. When passing the critical line N cr i t (£l) for the 
phase transition, the moment of inertia starts decreas- 
ing and the internal torque term K- mt changes sign. This 
leads to a narrow dip for Ct(N) in the vicinity of this 
line. As a result, the phase diagram gets overpopulated 
for N ~ A^ cr it(Sl) and depopulated for N ~ N cl -i t (fl) up 
to the second maximum of I(N, f2) close to the black- 
hole line N max (Q). This population gap marks the region 
of hybrid stars in the phase diagram and is a measur- 
able. Moreover, as we have seen population clustering 
of compact stars at the phase transition line could be a 
signal for the occurrence of stars with high density mat- 
ter cores and a measure for obtaining constraints on the 
EoS at high densities. In contrast to this scenario, in the 
case without a phase transition, the moment of inertia 
could at best saturate before the transition to the black 
hole region and consequently would also saturate. This 
would entail a smooth population of the phase diagram 
without a pronounced structure |u) . 

As a strategy for the investigation of high density phase 
transitions in compact stars we suggest to perform a sys- 
tematic observation of LMXBs for which the discovery 
of strong and remarkably coherent high-frequency QPOs 
with the Rossi X-ray Timing Explorer has provided new 
information about the masses and rotation frequencies 
of the central compact object |Tl]Jl^ , |30| |. If, e.g., the re- 
cently discussed period clustering for Atoll- and Z-sources 
P.p9[ will correspond to objects in a narrow region of 
masses well below the maximum mass limit, this could 
be interpreted as a signal for the high density phase tran- 
sition to be associated with a fragment of the critical line 
in the phase diagram for rotating compact stars flC|| . 

In the present work we have shown that the exis- 
tence and the shape of the suggested population gap for 
LMXBs in the phase diagram will signal the occurence 
of a phase transition in the QCD EoS and constrains its 
properties at high densities. 
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